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(54) Title: A. CODE SYNCHRONIZATION UNIT AND METHOD 

(57) Abstract 

A pilot acquisition unit for code division multiple 
access (CDMA) communication systems provided which 
includes a fast Hadamard transform (FHT) unit and a 
Sadamard processing unit. The FHT unit ^rrmnc* 
mTqualLty, in 'accordance with a mra of each of a set 
of possible pseudo-random number <PN) loadmgs and the 
ore-Hadamard processing unit generates a vector u per set 
of PN loadmgs. The vector u derines a quahty metric of a 
Scd rXt signal with the set of possible FN loadings, 
the pie-Hadamard processing unit providing the vector u to 
the FHT unit. 
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A CODE SYNCHRONIZATION UNIT AND METHOD 
FIELD OF THE INVENTION 

The present invention relates to digital communication systems 
generally and to a method and apparatus for acquisition of digital communication 
signals. 

BACKGROUND OF THE INVENTION 

Digital communication systems transmit and receive signals which 
have digital information therein. Typically, such signals include the data to be 
transmitted plus additional portions needed to ensure accurate communication, 
such as synchronization signals (to synchronize the receiver with the transmitter) 
and error correcting codes (to ensure that the received data has not been 
corrupted and to correct at least part of any corrupted data). 

There" are many types of digital communication systems. One 
5 common one is that of a spread spectrum system. A conventional direct 
sequence spread spectrum signal can be viewed as the result of mixing a 
narrowband information-bearing signal with an informationless wideband (and 
constant envelope) "spreading" signal. If B, and B p denote the bandwidths of the 
information-bearing signal and the spreading signal, respectively, then the 
20 "processing gain" available to the receiver is G = B/B p . The receiver synchronizes 
the incoming signal to a locally generated version of the spreading signal and 
mixes the received signal with the locally generated spreading signal, thereby, 
removing the spreading signal from the received signal and "collapsing" the signal 
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to the "information bandwidth" B,. 

The spread^ sianai is a coded sequence of some «. such 
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Switches 16, also known as taps, have predetermined states h k and 
are either closed (h k = 1) or open {\ = 0). The initial switch h 0 is always closed 
and provides the output c t directly to the m-1th flip-flop 12. The pseudo-random 
code sequence p[t] is composed of PN symbols, the duration of each of which is 
termed a "chip". Each symbol of the sequence is defined by: 
Equation 1 

In order to synchronize the local version of the spreading signal with 
the original version, the transmitting unit additionally transmits the pilot signal, 
containing the code sequence. For simplicity, we assume that the transmitted 
signal is binary phase shift keying (BPSK) modulated. 

The local unit then synchronizes its local code generator to the pilot 
signal after which the local unit can despread the received information bearing 
signals. The pilot signal is also utilized to track variations in the transmission 
5 channel. 

The received signal, after being down converted to a baseband signal 
and sampled at the output of a matched filter at a rate of one sample per chip, is 
denoted by t =■- ,-2,-WA- . The received signal consists of the pilot 

signal and the user data signals, both of which are transmitted by the transmitting 
20 unit, and interference terms caused by thermal noise and by signals transmitted 
by adjacent transmitting units. 

For the purpose of acquiring the initial synchronization, only the pilot - 
signal P ilot[t] is of interest. For a BPSK signal, the pilot signal may be 
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represented by: 
Equation 2 

r-i 

where p°[t] is a PN sequence, a^* 1 is the channel gain of the / -th 
signal reflection (called a "finger"). F denotes the number of fingers and v n 
denotes the residual frequency drift after baseband down-conversion. Now, 
consider only the most significant finger (the one with the largest a,) and denote it 

by ae J * . Also, denote the contribution of all other fingers, the user data signals 

and other interferences by n[t] . Then R[t] is represented by: 

Equation 3 

R[t]=ap°[t¥ M +n[t] 

The acquisition problem is how to efficiently obtain the phase of the PN 
sequence (i.e. the current loading of the PN generator 10) given some 
measurement record R[t] t = L2,-~,N. 
5 Solutions to the acquisition problem are described in the book CDMA: 

Principles of Spread Spectrum Communication , by A. J. Viterbi, Addison-Wesley, 
1995, in particular in section 3.4.3, pp. 58 - 59. The book is incorporated herein 
by reference. 

The direct approach is. to enumerate over all possible 2" - 1 phases of 
20 the PN sequence (there are 2 M possible initial loadings, but the zero loading is 
illegal since it produces an all zero sequence) and select the one which is optimal, 
with respect to some criterion. This approach is computationally and time 
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intensive due to the large number of possible PN loadings. 

A possible refinement of this approach, discussed in the book CDMA: 
Principles of Spread Spectrum Communication , is to obtain the phase by using a 
two-stage (dual-dwell) search procedure, where the first stage enumerates over 
5 all possible PN phases and passes only those phases with metric values that are 
above some pre-specified threshold to the second stage. In the second stage, 
each phase hypothesis is examined more thoroughly (i.e., using a more 
computationally intensive criterion) in order to decide whether it is the true PN 
phase or not. The dual-dwell procedure is faster than the direct approach but still 
10 takes a significant amount of time. 

When there is no frequency drift in the received samples R[t] (i.e., 
oj 0 =0), the optimal metric, under a white Gaussian noise assumption, whose 
absolute value needs to be maximized is a Maximum Likelihood metric, as 
follows: 

15 Equation 4 

metric = £ 

where ltfl],il[2], "« the block of sampled data, sampled at the 

rate of one sample per chip, and p[t] is one possible PN sequence. In the dual 
dwell procedure, the size N of the block is relatively small in the first phase and 
20 larger in the second phase. 

If the data might have a frequency drift, the metric should be 
insensitive to frequency drifts. The following differential metric has been 
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suggested by M. H. Zarrabizadeh and E. S. Souza in the article "Analysis of a 
Differentially Coherent DS-SS Parallel Acquisition Receiver", IEEE Proceedings of 
the 45 th Vehicular Technology Conference , Vol. 2, pp. 2/1 - 275, 1995 (the article 
is incorporated herein by reference): 

5 Equation 5 

metric = Y,z[l]z{l -l] 

where 

z[l] = Y,R[lN c +t}p[lN c +t} 

N c is the number of chips used for the coherent summation (e.g. the 
10 number of chips per symbol which is 64 in the IS-95 CDMA standard), and N s is 
the number of z[ ] variables used for creating the final metric. For example, N/\s 
small (e.g. 5) for the first phase (dwell) and larger (e.g. 10) for the second phase. 

The fallowing articles and patents discuss transform-domain methods 
for soft decoding of PN loadings and error correcting codes in general when 
15 BPSK signaling is used. The articles are incorporated herein by reference. 

V. V. Lasev and V. D. Dvornikov, "Determination of the 
Phase of a Pseudorandom Sequence From its Segment Using Fast 
Transforms" , Radio Engineering and Electronic Phvsics , Vol. 26, No. 8, 
pp. 61 - 66 T Aug. 1981; 

20 M. Cohn and A. Lempel, "On Fast M-Sequence Transforms'^ . 

IEEE Transactions on Information Theory , pp. 135-137, 1977; 
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V. V. Losev and V. D. Dvornikov, "Recognition of Address 
Sequences Using Fast Transformations", Radio Fnaineerinq and 
Electronic Physics , Vol. 28, No. 8, pp. 62-69, Aug. 1983; 

S. Z. Budisin, "Fast PN Sequence Correlation by Using 
FWT", IEEE Proceedings of the Mediter ranean Electrotechnical 
Conference (MELECONV Lisbon, Portugal, April 1989, pp. 513-51.5; 

Y. Be'ery and J. Snyders, "Optimal Soft Decision Block 
Decoders Based on Fast Hadamard Transform" , IFFE Transactions on 
information Theory , Vol. 32, 1986, pp. 355-364; and 
U.S. Patent 5 ,463,657 to Rice. 
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SUMMARY OF THE PRESENT INVENTION 

It is an object of the present invention to provide a novel and relatively 
fast method and apparatus for synchronization to a pilot signal, particularly for 
CDMA systems. 

It is a further object of the present invention to provide a method and 
apparatus for synchronization to a pilot signal which has a frequency drift therein. 

It is a still further object of the present invention to provide a method 
and apparatus for soft decoding an error correcting code when frequency drift is 
present. 

There is therefore provided, in accordance with a preferred 
embodiment of the present invention, a pilot acquisition unit for code division 
multiple access (CDMA) communication systems which includes a fast Hadamard 
transform (FHT) unit and a pre-Hadamard processing unit The FHT unit 
determines the quality, in accordance with a metric, of each of a set of possible 
pseudo-random number (PN) loadings and the pre-Hadamard processing unit 
generates a vector u per set of PN loadings. The vector u defines a quality metric 
of a received pilot signal with the set of possible PN loadings, the pre-Hadamard 
processing unit providing the vector u to the FHT unit. 

Moreover, in accordance with a preferred embodiment of tine present 
3 invention, the unit includes a partial possible PN loading generator for generating 
a series of partial possible PN loadings s E , wherein each partial possible PN 
loading s 5 defines one set of possible PN loadings. 
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Further, in accordance with a preferred embodiment of the present 
invention, the unit includes a dual dwell unit for selecting the PN loadings having 
the metric values above a predetermined threshold from among the PN loadings 
selected by the local PN loading selector, for determining a second metric for 
each of the selected PN loadings and for selecting the PN loading from among 
the selected PN loadings with the best value for the second metric. 

Additionally, in accordance with a preferred embodiment of the present 
invention, the pre-Hadamard processing unit comprises a local PN generator and 
a u vector generator which performs the following steps: 

loads a local PN generator with an initial PN loading; 

loops on each of the datapoints of the received pilot signal and per 

loop: 

combines one partial possible loading s E with a datapoint of the 
received pilot signal and with a PN loading produced by the local PN generator 
5 thereby to update the u vector, 

steps the local PN generator to produce another PN loading; and 
provides the resultant u vector to the FHT unit 

Alternatively, for received signals with frequency drift therein, the 
pre-Hadamard processing unit comprises a local PN generator and a u vector 
20 generator which performs a similar set of steps as follows: 

loops over plurality of drift loop values, the step of looping including the 
steps of a) loading a local PN generator with a different initial PN loading per loop 
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value and b) generating an input signal which is insensitive to drift from the 
received pilot signal; 

loops on each of the datapoints of the input signal, .the second step of 
looping including the steps of a) combining one partial possible loading s £ with a 
datapoint of the input signal and with a PN loading produced by the local PN 
generator thereby to update the u vector and b) stepping the local PN generator 
to produce another PN loading. 

The remaining steps are the same as for the non-frequency drift case. 
Still further, in accordance with a preferred embodiment of the present 
invention, the unit can include a dual dwell unit which performs a further metric 
calculation on those PN loadings which produce a metric above a predefined 
. threshold. 

The present invention is operative for all digital communication 
systems (not just CDMA) which have frequency drifts therein and can also be 
5 implemented, as described and claimed hereinbelow, for signals encoded with 
error correcting codes. For the latter, the local PN generator is replaced with a 
generating matrix. 
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BRIEF DESCRIPTION OF THE DRAWINGS AND APPENDICES 

The present invention will be understood and appreciated more fully 
from the following detailed description taken in conjunction with the drawings and 
appendices in which: 

Fig. 1 is schematic illustration of a prior art pseudo-random number 
(PN) generator, 

Fig. 2 is a block diagram illustration of a pilot acquisition unit, 
constructed and operative in accordance with a preferred embodiment of the 
present invention; 

Fig. 3 is a schematic illustration of a pre-Hadamard processor forming 
part of the pilot acquisition unit of Fig. 2; 

Fig. 4 is a flow chart illustration of a method of operating the 
pre-Hadamard processor of Fig. 3 for signals with no frequency drift; 

Fig. 5 is a flow chart illustration of a method of operating the 
s pre-Hadamard processor of Fig. 3 for signals with frequency drift; 

Fig/ 6 is a block diagram illustration of a decoder for data encoding with 
error correction codes, constructed and operative in accordance with a further 
preferred embodiment of the present invention; 

Fig. 7 is a schematic illustration of a pre-Hadamard processor forming 
20 part of the unit of Fig. 6; 

Fig. 8 is a flow chart illustration of a method of operating the 
pre-Hadamard processor of Fig. 7; 
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Appendix A provides the mathematical basis for the pre-Hadamard 
processor of Figs. 2, 3 and 4 operating on quaternary phase shift keying (QPSK) 
signals having no frequency drift; 

Appendix B provides the mathematical basts for the pre-Hadamard 
processor of Fig. 5 operating on QPSK signals having an unknown frequency 
drift. 

Appendix C provides the mathematical basis for the pre-Hadamard 
processor of Figs. 2, 3 and 4 operating on binan/ phase shift keying (BPSK)- 
signals having no frequency drift; and 

Appendix. D provides the mathematical basis for the pre-Hadamard 
processor of Fig. 5 operating on BPSK signals having an unknown frequency drift. 
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DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

The pilot acquisition unit of the present invention considers each 
possible pseudo random number (PN) loading s (i.e. set of flip-flop values a,) and 
determines the value of a metric, metric[s] for it! The present invention then 
> reviews the set of metrics[s] and selects the PN loading which is associated with 
the "best" (e.g. largest in absolute value) metricfs]. The selected PN loading is 
the detected current PN loading of a PN generator on the transmitting unit with 
which the data was encoded. 

For CDMA systems which have quaternary phase shift keying (QPSK) 
o signals, the transmitted signals are complex and thus, the PN sequence is also 
complex: p*[t] = + The complex PN sequence is generated by two 

PN generators, one for the in-phase sequence p ; V]and one for the quaternary 
sequence p Q Q [t] . For such CDMA systems, the metric for obtaining p*[t] might 
be: 

15 Equation 6 

n * 
metritis] = 2>['lP/U] = Z *M(~ ^ 

/-I /ol 

where c t is the output of the in-phase PN generator and p,[t] is the 
in-phase portion of the complex QPSK PN sequence and is a ftinction of the 
in-phase PN loading s. As discussed hereinbelow, p^lt] is obtained from the 
20 estimated in-phase PN loading s. 

It will be appreciated that metric[s] of Equation 6 given above is valid* 
only if there is no frequency drift in the received data R[t]. 

13 
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It will further be appreciated that while the present invention is 

described with respect to CDMA QPSK signals, it is also operative with respect to 

other digital communication systems. 

The values of row vector metrics[s] (one element for each possible 

value of s) are generated via a fast Hadamard transform (FHT) by noting that 

metrics[s] can be written as follows: 

Equation 7 

metrics[s] = u • H m 

where H m is the Hadamard matrix and u is the input vector to the fast 
o Hadamard transform. The association of metrics(s] with the Hadamard matrix 
and the construction of the input vector u from the received samples are derived 
in detail in Appendix A for a CDMA, QPSK input pilot signal. Fast Hadamard 
transforms are discussed in the book Fast Transforms, Algorithms, Analysis, 
Applications , by D. F. Elliot and K.R. Rao, Academic Press, New York, 1982. The 
15 book is incorporated herein by reference. 

Reference is now made to Fig. 2 which illustrates, in block diagram 
format, the pilot acquisition unit of the present invention, to Fig. 3 which details a 
pre-Hadamard processor useful in the pilot acquisition unit of Fig. 2 and to Fig. 4 
which illustrates, in flow-chart format, the method of operating the pre-Hadamard 
20 processor of Fig. 3. 

The acquisition unit comprises a fast Hadamard transform (FHT) unit 
20, a pre-Hadamard processor 22, a partial, possible PN loading s E generator 24,. 
a local PN loading selector 26" and a global PN loading selector 28. Partial, 
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possible FN loading s £ generator 24 is typically a counter which provides the 
count value as the partial PN loading s E 

As described in more detail hereinbelow, pre-Hadamard processor 22 
produces the Hadamard input vector u for all of the loadings s which have the 
current partial, possible loading s E in common, given the received data R[t] and an 
initial loading h of a PN generator 19 (Fig. 3) forming part of pre-Hadamard 
processor 22. Initial loading h, defined as (h^. ... h„). is produced by providing 
the values of the taps 16 into their corresponding flip-flops 12, where a Q receives 
the value of h^ etc. 

FHT unit 20 performs a fast Hadamard transform on the Hadamard 
input signal u and produces therefrom the vector metrics[s] for all of the loadings 
s which have the current partial, possible loading s E in common. Local PN 
loading selector 26 selects the PN loading s, associated with the maximal 
component of |metrics[s]|. The process is repeated for all partial, possible 
5 loadings s s and global PN loading selector 28 selects the detected PN loading s 
from among those loadings s, produced by local PN loading selector with the 
largest value of |metric[s]j. 

IS-95 CDMA systems have two local PN generators, an in-phase PN 
generator and a quaternary PN generator, which are tied to each other and. the 
20 in-phase PN generator influences the sequence of the quaternary PN generator. 
The opposite is not true.. To synchronize both local PN generators with the 
transmitting PN generators, the local PN generators are initialized with their initial - 
loadings and are stepped together until the in-phase PN generator achieves the 

15 
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selected loading s. The quaternary PN generator will have achieved its 
appropriate loading. 

As shown in Fig. 3, pre-Hadamard processor 22 comprises a local 
pseudo-random number generator 19, similar to PN generator 10 of Fig. 1, a 
Hadamard vector u register 30, a summer 32, a scalar multiplier 34 and a 
XOR-AND unit 36. 

In accordance with the article by Be'ery and Snyders discussed herein 
in the Background and as detailed in Appendix A, PN generator 19 is divided into 
two sections, an external section E of length m-Q incorporating the m-Q flip-flops 
12 having values a Q to and an internal section I of length Q incorporating the 
Q flip-flops 12 having values a^ to a^. The internal section creates an internal 
'vector g* t and the external section creates an external vector gf where the vectors 
g are defined by: 



and the a! are the values in the flip-flops 1 2 at time t 
It will be appreciated that the vector g, = (g,',gf ) is one possible state 
of the local PN generator 19 while PN loading s is the PN loading with which the 



20 output c t of the PN generator 19 is a function of the PN loadings s and g t as 



Equation 8 




5 



received pilot signal was generated. Furthermore, it will be appreciated that the 



follows: 



Equation 9 
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where < > denotes a XOR-AND operation and where XOR is denoted 

by ©: 

<x,y>=(;c 0 AND y 0 )©( x i ^) e ---®<X-« 
5 The value of internal vector g[ defines an address within Hadamard 

vector u register 30, where register 30 contains 2 Q memory ceils. Arrow 40 which 
points from internal vector g| to the address it defines, labeled 42. 
Pre-Hadamard processor 22 removes the value stored, in address 42 and 
provides the value to summer 32. 
10 The external vector gf is utilized, in combination with the possible, 

partial PN loading s E , to determine the sign of the datapoint R[t]. It is noted that a) 
the partial PN loading s E is of the same length m-Q as the external vector gf and 
b) the external vector gf and the partial PN loading s E are binary vectors of 1's 
and O's. Specifically, the operation performed is: 

15. Equation 10 

sign = (-1) 

Scalar multiplier 34 multiplies the datapoint R[t] by the value of sign 
and the result is added, in summer 32 f to the component of the Hadamard input 
vector u removed from address 42. The output of summer 32 is then inserted 
20 , back into address 42. 

As indicated in Fig. 4, pre-Hadamard processor 22 repeats the 
above-described operations for each of the N values of R[t]. Initially (steps 50 and 

17 
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51), pre-Hadamard processor 22 zeros the Hadamard vector u register 30 and 
loads the PN generator 19 with its initial state vector h. In step 52 T pre-Hadamard 
processor 22 loops over the N values of R[t] where, for each value of R[t], the new 
value for the relevant component of u is determined (step 54) after which the PN 
generator 19 is stepped (step 56) to produce new values for internal vector g ! t 
and external vector gf . 

Once loop 52 has completed, the u vector has been produced and, 
therefore, can be sent (step 58) to FHT unit 20 for determining the values of 
metrics[s]. The process begins again at step 50 by resetting the Hadamard 
vector u and PN generator 19 to their initial states. 

It will be appreciated that FHT unit 20 operates once per partial PN 
loading s E while pre-Hadamard processor 22 repeats its operations N times per 
partial PN loading s E . The number of partial PN loadings s £ is 2"^° where Q is 
chosen to balance between the number of operations performed by 
5 pre-Hadamard processor 22 and the number of operations performed by FHT unit 
.20. 

It will further be appreciated that the fast Hadamard transform 
performed by FHT unit 20 performs a series of addition operations only. As a 
result, the pilot acquisition unit of the present invention performs addition 
20 operations only (there are no real number multiplications since XOR-AND unit 36 
performs only XOR-AND operations and scalar multiplier 34 only produces a sign 
change). Since the number of addition operations is relatively low, the pilot' 
acquisition unit of the present invention performs the pilot synchronization 

18 
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operation faster than in the prior art. Just how much faster depends on the 
selection of Q. 

For example, there might be N=840 samples in the received signal 
R[t], the length m of the PN generators might be 15, and the split value Q might 
be 12. 

The relatively fast acquisition is particularly useful for CDMA systems, 
such as for cellular telephony, where initial synchronization needs to be acquired 
as quickly as possible. 

It will be appreciated that the equation for metric[s] given above in 
Equation 6 is valid only if there is no frequency drift in the received data R[t]. 
However, this is rarely the case. As discussed in Appendix B, the baseband 
down-conversion process is not ideal and some residual frequency drift will 
always be present (typically due to clock rate mismatch between the transmitter 
and receiver). 

5 One possible solution is to enumerate over all possible frequency drifts 

and to eliminate the drift from the received data for each such hypothesis. Then 
the method that was presented above may be applied on the transformed data. 
Another possibility is to utilize a metric which is insensitive to the frequency drift 
One such metric is a multi-differential metric, provided hereinbelow for QPSK 

20 signals (the metric for BPSK signals is provided in Appendix D): 

Equation 11 

m e mc[s]=tl^[fK['l 
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where 



and the number L of differentials is a small, predetermined number. 
For simplicity, L is less than or equal to m but other values are possible. 

5 As derived in Appendix B, the multi-differential metric is determined in 

a manner similar to that of the non-frequency drift metric and thus, the pilot 
acquisition unit of Fig. 2 can be utilized to determine the PN loading associated 
with the best multi-differential metric. However, for this, multi-difFerential 
embodiment, the pre-Hadamard processor 22 of Fig. 2 is operated in accordance 

10 with Fig. 5 T to which reference is now made. 



processor 22 per partial PN loading s E . Initially (step 70), pre-Hadamard 
processor 22 zeros the Hadamard vector u register 30. Following the preparation 
of register 30, pre-Hadamard processor 22 begins a loop 72 over the possible 
15 values of / . For each value of / f pre-Hadamard processor 22 generates (step 
74) the Mh input data loading R, as per Equation 11.. Pre-Hadamard processor 
22 also generates (step 76) the Mh loading h i of PN generator 19 and loads it 
(step 78) into PN generator 19; The / -th loading h, is defined by: 



Similar to Fig. 4, Fig. 5 illustrates the operations of pre-Hadamard 
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Equation 12 
h, = h®z, 
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where h is the initial loading of PN generator 19 and © indicates the 
XOR operation. 

Pre-Hadamard processor 22 then determines the values of u register 
30 in loop 80 as in the previous embodiment with the following exceptions: 

i) PN generator 19 is loaded with its / -th loading h, rather than h; and 

ii) the input data sequence R, replaces the previous input data 
sequence R. 

Specifically, u register 30 is updated (step 82) after which the PN 
generator 19 is stepped. 

After both loops 72 and 80 are finished, the vector u stored in register 
30 is provided (step 86) to FHT unit 20. The remaining operations of the pilot 
acquisition unit are as before. Pre-Hadamard processor 22 repeats its operations 
per partial loading s E , local selector 26 selects the best solution s, per partial 
loading s E , and global selector 28 selects the best overall loading s. 

For example, there might be N=2560 samples in the received signal 
R[t], the length m of the PN generators might be 15, the split value Q might be 12 
and the number of differentials L might be 8. 

It will be appreciated that the pilot acquisition unit of the present 
invention can be utilized for systems with and without frequency drift, the only 
difference being the initial loading of the PN generator and the differential or 
non-differential input data. It is noted that all digital communication systems using 
PN generated codes, such as CDMA systems, other spread spectrum systems 
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and systems which' add error correcting codes to transmitted data, regardless of 
the format (QPSK t BPSK, etc.) of the transmitted data, typically have some 
frequency drift therein and thus, the second embodiment of Fig. 5 is typically 
applicable. 

The pilot acquisition unit of the present invention can be operated 
within a dual dwell scheme, as follows. At the first stage (dwell), the unit of Fig. 2 
detects all PN loadings which produce a metric value (with or without frequency 
drift) above some pre-selected threshold. At the second dwell, a further unit (not 
shown) calculates the prior art metric defined either by Equation 4 (no frequency 
drift) or by Equation 5 (with frequency drift), where, in both, the signal p[t] is 
replaced by the signal p*[t] in order to deal with a QPSK signal. 

The threshold is determined as follows. Let the a-posterio'ri variance of 
the metric (given the data measurements) under a random loading be denoted by 
cr (each metric has its own value of cr). At each dwell, only the hypotheses that 
are above t-c (i.e. |mem*c(s)| > t - cr ) are passed on to the next stage (either the 
hypothesis is passed on to the second dwelt or it is selected as a successful 
synchronization). For example, in a frequency drift situation, t might be set as 
t = 4.0 for the first dwell (using the metric of Equation 1 1 ) and t = 55 for the 
second dwell (using the metric of Equation 5). 

Note that the a-posteriori variances of the various metrics are given 

by: 
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2 



2j:\\R[t}i 



for the metric defined by Equation 4 



4l IWt-iytr.+ti 



for the metric defined by Equation 5 




for the metric defined by Equation .6 



L N 



a' 



for the metric defined by Equation 11 



As mentioned hereinabove, the IS-95 specification for CDMA defines a 
complex PN loading, p°[f] = P/ a [*] + y>g[f], which is generated by two PN 
generators. The first PN generator generates p° r [t]. The second PN generator 
generates Both PN generators are initialized at the beginning of the 

transmission. 

In order to increase the periods of the PN sequences from 2 15 - 1 to 
2 15 , and to balance the number of Q's and Vs in these sequences, the following 
non-linear mechanism is employed. Whenever the pattern 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,1 is detected in the first sequence, an additional 0 bit is 
inserted into the output of both sequences. However, the present invention does 
not collapse in the presence of this non-linearity but, instead, is degraded slightly. 

It will be appreciated that the present invention is not limited by what 
has been described hereinabove and that numerous modifications, all of which fall 
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within the scope of the present invention, exist. For example, while the present 
invention has been described with respect to CDMA systems, it can be 
implemented in other digital communication systems. In particular, the present 
invention incorporates all implementations of code synchronization in the 
presence of frequency drifts, whether in conjunction with a pilot signal or not 

Reference is now made to Figs. 6, 7 and 8 which illustrate a decoder 
for messages encoded with error correcting codes which utilizes the concepts of 
the present invention. For this embodiment, let s be the message and let p[t] be 
the encoded version of the message s (of length m) which is the sequence to be 
transmitted. Once again, the received and sampled signal is R[t]. 

To encode the message s, a generating matrix G, with columns g tt is 
utilized, where: 

c?-[«f «r -] 

A value c, is a function of the message s and the t-th generating vector 
g, (of length m) and the sequence p[t] is produced from the value c, as follows: 

e, =< s,g, > 

pM=Hr=(-ir ,,> 

Thus, the sequence p[t] to be transmitted has the same structure as in 
the previous embodiments (see Equations 6 and 9) although it is formed from 
different components. 

The metric to be maximized must be insensitive to frequency drift". For 
BPSK modulation, the metric is: 
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metric&s) = Y £,[*]y[r] 
/-i 

where, as in Appendix D, 

As in Appendix D, we assume l^m and we have 

«,M = « 2 M P°MO>°[' - '])" a**) + *1 « M^'M' " 'fl") + *] 
where 77M denotes the contribution ; of the noise terms. Now, for BPSK 
modulation, 

and 
Thus: 

p o [^°^-n=(-ir s ' e5 *^ 

Fig. 6 illustrates the decoder of the present invention for data encoded 
with error correcting codes. It has a similar structure to that of the pilot acquisition 
5 unit and thus, similar elements carry similar reference numerals. 

The decoder comprises FHT unit 20, a pre-Hadamard processor, 
labeled 90, partial, possible s E generator 24, a local message selector 100 
operating similarly to local PN loading selector 26 and a global message, selector 
102 operating similar to global PN loading selector 28. 
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Similar to that of the pilot acquisition unit of the previous embodiments, 
pre-Hadamard processor 90 produces the Hadamard input vector u for all of the 
messages s which have the current partial, possible message s E in common, 
given the received data R[t], 
5 FHT unit 20 performs a fast Hadamard transform on the Hadamard 

input signal u and produces therefrom the vector metrics[s] for ail of the . 
messages s which have the current partial, possible message s E in common. 
Local message selector 100 selects the message s, associated with the maximal 
component of |metrics[s]|. The process is repeated for all partial, possible 
10 message s E and global message selector 102 selects the detected message s 
from among those messages s, produced by local message selector with the 
largest value of |metric[s]l. 

Fig. 7 illustrates the pre-Hadamard processor 90 which is similar to " 
pre-Hadamard processor 22 of Fig. 3 in that it comprises Hadamard vector u 
15 register 30, summer 32, scalar multiplier 34 and XOR-AND unit 36. However, 
pre-Hadamard processor 90 comprises storage unit 104, storing the generating 
vectors of generating matrix G and a XOR unit 106 instead of the local PN 
generator 19 of pre-Hadamard processor 22. Storage unit 104 and XOR unit 106 
together produce the internal and external vectors, labeled g[ t/ and gf[ { 
20 respectively, which the Hadamard vector u register 30 and XOR-AND unit 36 
require. 

Specifically, XOR unit 106 generates a combination generating vector" 
g u from two vectors g t and g w which are stored in the storage unit 104. The latter 
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vector is / vectors away from the former where / is as defined hereinbelaw. 

As in the previous embodiments, the combination generating vector g u 

is divided into internal and external vectors, g ! t t and g* , where the internal vector 
g' t , contains Q components of the combination vector g u and the external vector 
gf, contains m-Q components of the combination vector g M . XOR-AND unit 36 



combines the external vector g* with the partial possible message s E , as 
described hereinabove in equation 10, and the multiplier 34 combines the result 
with the shifted received data R f [t] , defined hereinabove, 

As in the previous embodiments, internal vector g,' t/ is utilized to define 
an address within register 30. This is indicated by arrow 40 which points to the 
address, labeled 42. Pre-Hadarnard processor 90 removes the value stored in 
address 42 and provides the value to summer 32. 

Fig. 8 illustrates the operations performed by pre-Hadamard processor 
90. They are similar to those shown in Fig. 5 except that the operations on a pN 
5 generator are replaced with those on the generating matrix G. Specifically, 
pre-Hadamard processor 90 begins by zeroing (step 70) the Hadamard vector u 
register 30. Following the preparation of register 30, pre-Hadamard processor 90 
begins loop 72 over the possible values of I . For each value of I , pre-Hadamard 
processor 90 generates (step 74) the /-th input data loading R, and generates 
20 (step 1 10) the combined generation vector g u using the current value of / . 

Given combined generation vector g u , pre-Hadamard processor 90 
then determines the values of u register 30 in loop 80 as discussed hereinabove 

27 
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and update u register 30. 

After both loops 72 and 80 are finished, the vector u stored in register 
30 is provided (step 86) to FHT unit 20. 

it will be appreciated by persons skilled in the art that the present 
5 invention is not limited to what has been particularly shown and described 
hereinabove. Rather the scope of the present invention is defined only by the 
claims which follow after the Appendices: 
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Appendix A : ** 
Suppose that ai! transmitted signals are QPSK modulated, and let the 
received complex CDMA signal after baseband down-conversion, matched 
filtering and sampling be denoted by 4*] t =■•• -2 -1,0,U>'"« R [ £ ] c° nsists of 
5 the following components: 

1 .Pilot signal that is transmitted by the base-station. 
2.User data signals that are transmitted by the base-station. 
3. Interference terms including thermal noise and signals transmitted by 
adjacent base-stations. 
10 For the purpose of acquiring initial synchronization we are only 

interested in the pilot signal, pibt[t], which may be represented by 

/-i 

where p°[t] is a complex PN sequence composed of an in-phase PN 
sequence p*[t] and a quaternary PN sequence (j^M, p Q [t] = p?[t] + jpl[t])* 
15 a^'is the complex channel gain of the Mh finger, F denotes the number of 
- fingers., and a 0 denotes the residual frequency drift after baseband 
down-conversion. Now, consider only the most significant finger (the one with the 
largest a t \ and denote the contribution of all other fingers, and of the user data 
signals (component 2 above) and other interferences (component 3 above) by 
20 n[t]. Then R[t] is represented by: 
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Equation 13 

We assume that n[t] is a zero mean white noise term with variance 



The problem is how to efficiently obtain the phase of both PN 
sequences (i.e. the current loading of the PN generators) given some 
measurement record R[t] t = L2,- ,N • With no frequency drift, the metric for 
obtaining the in-phase PN loading is: . 
Equation 14 
n 

metric = J^R[t]p r [t] 

ml 

To calculate Equation 14 over all possible PN sequences p r [t], 
efficiently, we use a block-code soft decoding method, as follows hereinbelow. Let 
c, be the output bit produced by the PN generator or linear feedback shift (LFSR) 
5 at time t . Let the LFSR (Fibonacci form) have m cells and transition matrix given 
by 



M 



0 1 



A. 0 
1 0 



0 1 
0 0 



so 



that a (,) = a ( '-°A/ , where a"> -(af^-o^) is the state of the shift 
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register at time t . 

Now let the state of the LFSR at time / = lbe denoted by 
c-k s s ). Given the data measurements H[l],ii[2],-,ii[^] , we need to 
obtain s . Note that 

where c, =sM'M T = sgf »(s,8,)' . g,-h(M') r and h = (^,A„,. 1 .-a) 
(/lo s 1)t Note that g, is the state of the Galois form LFSR after t clocks, when 
initialized by h. The Galois form LFSR has transition matrix M T , i.e. 
a ( ' ) = a ( '- |> M r . 

10 Now, given a vector y = (y 0 .3V ^-i) • we define 

Hence, y = 2' 

Let A, = {t : &(g,) = ;) and = Z*M • Hence, 
me mc(s)=E(-l) (, - i " 0 ' )) "; 

y-o 

15 so that 

(m«^6- , (0)).-".«eO*(l'" , ( 2 "- l )))" u - i7 -- 
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where H,„ is the Hadamard matrix, defined by {H M ) tJ = ^ ^ ■ 
A crucial point in the algorithm is that H m may be calculated efficiently using the 
following recursion, 

5 Hence, for any vector n of dimension 2'" , if u is partitioned to two 

sub-vectors n, , u 2 of dimension 2~' each, i.e. n = (u„u 2 ) , we have 
Equation 15 

u • H„ + u 2 /f m _, , u,H,„., - u,^,) 

That is to say, the Hadamard transform (HT) of u may be obtained 
10 from the HT-s of u, , u, . Furthermore, Equation 15 can be recursively invoked on 
each of the smaller dimensional HT-s to produce the Fast Hadamard transform 
(FHT) algorithm. 

The system may be improved in terms of both computational time and 
memory requirements as follows. Let g, and s be partitioned as follows 

is S,=(g,'gf) ■-(»'■') 

where and a* have dimension O, and where gf and s E have 
dimension m-Q (1 denotes internal space, and E denotes external space). 
Hence, 

metrics) -t^-lt'^i-ir^m 
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Let, 

Aj = {t : b{°:)=j} j = 0,l,-2 Q -l 
and 

Equation 16 

4° = Z(-if' Bi!) *M 



10 



Then, 



7-0 



To calculate metric(s) efficiently we enumerate over all possible values 
of s s . For each value of s £ , we first calculate uf j = 0,1,- 2 Q -1 using 
Equation 16, and then apply the fast Hadamard transform. The estimated loading 
of the PN is that value of s = (V s E ) that maximizes |mernc(s)|j 2 . 
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1 Appendix B 

In most transmission systems, the received samples are subject to an 
unknown frequency drift (i.e., <u fl is non-zero). For this case, the suggested metric 
is 

5 Equation 17 

where 

For simplicity we assume I < m . Recall Equation 3, we have 

where 77M denotes the contribution of the noise terms. The 
approximation is due to the fact that the frequency drift is typically such that, for 
small I, (D 0 l «1- 
Now, 

In addition, 

where 

Equation 18 
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Let h, =h®z r Then p?[f]p?[r-/] = (-l) 8A ' lh,r . The metric is maximized 
for p l j[t] = (-l) ,w ' h ' . Hence, the algorithm that was presented in Appendix A may 
be applied to R t \t] , except that the initial loading of the PM is h, instead of k. We 
call this algorithm, the differential algorithm. We now suggest a multi-differential 
algorithm, which is an extension of the differential algorithm and is described in 
the text hereinabove with respect to Fig. 5. The suggested metric is 
Equation 19 

metrics) ^YtJ^RtUWlt] 

/-I r-l 

The metric is optimized for ^W = (-l) BA/<hr ■ Denote gf / = k,(M r )' • 

Then 

i«l /-I 

We proceed by using an algorithm, which is similar to the one that we 
used in the previous case (no frequency drift). Let g rJ and s be partitioned as 

follows 

where and s 7 have dimension O, and where g* and s £ have 
dimension m-O. (I denotes internal space, and E denotes external space). 
Hence, 
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10 



metri 



Let, 



15 



A^{f.b{ S ' l4 ) = j\ 7 = 0,1,-^-1 
£-&(■*) 

and 

Equation 20 
Then, 

^n C ( S )=S(-l) {, '^ ,y)> ^ 

To calculate metrics) efficiently we enumerate over all possible values 
of s E . For each value of s B , we first calculate uf j = 0,l,-2 Q -1 " sina . 
Equation 21, and then apply the fast Hadamard transform. The estimated loading 
of the PN is that value of s = (s' s E ) that maximizes \metric(sf . Basically, Fig. 3 
describes the algorithm, except that to create {««) the PN generator 10 needs to 
be reloaded and advanced L times {L <m), in order to create L sequences: 

for Z = 1,2,---,I. The fast Hadamatti transform routine needs to be 
applied only once, per each value of s* . The benefit of the multi-differential (over 
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the differential) system is that a smaller amount of data is required. 
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APPENDIX C 

BPSK modulation utilizes a singie binary sequence which results in a 
different metric. However, the derivation is similar, as will be discussed 
hereinbelow. 

Without frequency drift, the pilot signal has the form (provided in the 
Background as Equation 2): 
Equation 21 

As before, consider only the most significant finger and denote the 
contribution of all other fingers, the user data signals and other interferences by 
n[t] . Then R[t] is represented by Equation 3 (repeated here): 
Equation 22 

For BPSK modulation, p[t] and p°[t] are reals and not complex and 
5 thus, the metric is provided in Equation 4 (repeated here): 

Equation 23 

metric = X 4'1p(0 

The rest of the derivation is identical to that provided in Appendix A 
(after Equation 14) where, for this embodiment, p t [t] of Appendix A is replaced by 
20 pffl. 
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APPENDIX D 

For BPSK modulation with frequency drift, the metric is similar to that of 
Equation 19, as follows: 

metric/is) -Y.R,[t]p l [t] 

5 where, as in Appendix B, 

As in Appendix B, we assume I £ m and we have 
R t [t] = <r Re(/[r]0>°l< - § *»} + *l - ^ M^tA' " l #] + *l 
where i£t] denotes the contribution of the noise terms. Now, for BPSK 
10 modulation, 

te{p o {t}{p o [t-i\)'}=P 0 [t]p a [t-n 

and 

^3=(-ir hr P o [<-']=(-ir v 

where z, is defined by Equation 18. 
The metric is maximized for 



15 



The rest of the derivation is identical to that provided in Appendix A 
(after Equation 18) where, for this embodiment, pffl of Appendix B is replaced by 
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CLAIMS 

1 . A pilot acquisition unit for a code division multiple access (CDMA) 
communication systems, the pilot acquisition unit comprising: 

a. a fast Hadamard transform (FHT) unit for determining the 
5 quality, in accordance with a metric, of each of a set of possible 

pseudo-random number (PN) loadings; and 

b. a pre-Hadamard processing unit for generating a vector u per 
set of PN loadings, said vector u defining a quality metric of a 
received pilot signal with said set of possible PN loadings, said 

10 pre-Hadamard processing unit providing said vector u to said 

FHT unit. 

2. A pilot synchronization unit according to claim 1 and wherein said 
pre-Hadamard processing unit comprises a local PN generator. 

3. A pilot synchronization unit according to claim 2 and also comprising a 
15 partial possible PN loading generator for generating a series of partial 

possible PN loadings s E , wherein each partial possible PN loading s E 
defines one said set of possible PN loadings. 

4. A pilot synchronization unit according to any of the previous claims 
and also comprising a local PN loading selector for selecting the PN 

20 loading for which said FHT produces the best metric. 

5. A pilot synchronization unit according to claim 4 and also comprising a 
global PN loading selector for selecting the PN loading having the 
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best metric from among the PN loadings selected by said local PN 
loading selector. 

6. A pilot synchronization unit according to claim 4 and also comprising a 
dual dwell unit for selecting the PN loadings having the metric values 
above a predetermined threshold from among the PN loadings 
selected by said local PN loading selector, for determining a second 
metric for each of said selected PN loadings and for selecting the PN 
loading from among the selected PN loadings with the best value for 
said second metric. 
7. A pilot synchronization unit according to claim 3 and wherein said 

pre-Hadamard processing unit comprises u vector generating means 

including: 

a. initializing means for loading said local PN generator with an 
initial PN loading; and 

b. loop means operating on each of the datapoints of said 
received pilot signal, said loop means including; 

i) combination means for combining one partial possible 
loading s E with a datapoint of said received pilot signal 
and with a PN loading produced by said local PN 
generator thereby to update said u vector, and 
. ii) PN stepping means for stepping said local PN 
generator. 



41 



PCT/IL98/00067 

WO 98/4467(1 

8. A pilot synchronization unit according to claim 3 wherein said received 
pilot signal has a frequency drift therein and wherein said 
pre-Hadamard processing unit comprises u vector generating means 
including: 

5 a. drift loop means operating on a plurality of loop values for 

loading said local PN generator with a different initial PN 
loading per loop value and for generating an input signal from 
said received pilot signal which is insensitive to drift; and 
b. datapoint loop means operating on each of the datapoints of 
10 said input signal, said loop means including: 

i) combination means for combining one partial possible 
loading s E with a datapoint of said received pilot signal 
and with a PN loading produced by said local PN 
generator thereby to update said u vector, and 
15 ji) pn stepping means for stepping said local PN 

generator. 

9. A code synchronization unit for synchronizing to received encoded 
signals having frequency drift therein, the code synchronization unit 
comprising: 

a. a fast Hadamard transform (FHT) unit for determining the 
quality, in accordance with a metric, of each of a set of possible 
pseudo-random number (PN) loadings; and 



20 
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b. a pre-Hadamard processing unit for generating a vector u per 
set of possible PN loadings, said vector u defining a quality 
metric of said received encoded signal with said set of possible 
PN loadings, said pre-Hadamard processing unit providing said 
vector u to said FHT unit. 

10. A code synchronization unit according to claim 9 and wherein said 
pre-Hadamard processing unit comprises a local PN generator. 

1 1 . A code synchronization unit according to claim 10 and also comprising 
a partial possible PN loading generator for generating a series of 
partial possible PN loadings s E . wherein each partial possible PN 
loading s E defines one said set of possible PN loadings. 

12. A code synchronization unit according to any of claims 10 - 11 and 
also comprising a local PN loading selector for selecting the PN 
loading for which said FHT produces the best metric. 
5 1 3. A code synchronization unit according to claim 12 and also comprising 
a global PN loading selector for selecting the PN loading having the 
best metric from among the PN loadings selected by said local PN 
loading selector. 

14. A code synchronization unit according to claim 12 and also comprising 
a dual dwell unit for selecting the PN loadings having the metric 
values above a predetermined threshold from among the PN loadings 
selected by said local PN loading selector, for determining a second 
metric for each of said selected PN loadings and for selecting the PN 
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loading from among the selected PN loadings with the best value for 
said second metric. 
15. A code synchronization unit according to claim 11 wherein said 

pre-Hadamard processing unit comprises u vector generating means 

including: 

a. drift loop means operating on a plurality of loop values for 
loading said local PN generator with a different initial PN 
loading per loop value and for generating an input signal from 
said received pilot signal which is insensitive to drift; and 
b. datapoint loop means operating on each of the datapoints of 
said input signal, said loop means including: 

i) combination means for combining one partial possible 
loading s E with a datapoint of said received pilot signal 
and with a PN loading produced by said local PN 
generator thereby to update said u vector; and 
,i) PN stepping means for stepping said local PN 



15 

ii 



generator. 

16. An error correcting code decoder for decoding signals having 
frequency drift therein, the decoder comprising: 
20 a . a fast Hadamard transform (FHT) unit for determining the 

quality, in accordance with a metric, of each of a set of possible 
messages; and 
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b. a pre-Hadamard processing unit for generating a vector u per 
set of possible PN loadings, said vector u defining a quality 
metric of said received encoded signal with said set of possible 
messages, said pre-Hadamard processing unit providing said 
vector u to said FHT unit 
17. A decoder according to claim 16 and also comprising a partial 
possible message generator for generating a series of partial possible 
messages s E , wherein each partial possible message s £ defines one 
said set of possible messages. 

18. A decoder according to any of claims 16 - 17 and also comprising a 
local message selector for selecting the message for which said FHT 
produces the best metric. 

19. A decoder according to claim 18 and also comprising a global 
message selector for selecting the message having the best metric 

5 from among the messages selected by said local message selector. 

20. A decoder according to claim 18 and also comprising a dual dwell unit 
for selecting the messages having the metric values above a 
predetermined threshold from among the messages selected by said 
local message selector, for determining a second metric for each of 

20 said selected messages and for selecting the message from among 

the selected messages with the best value for said second metric. 
21. A decoder according to claim 16 wherein said pre-Hadamard 
processing unit comprises u vector generating means including: 
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a. drift loop means operating on a plurality of loop values I for 
creating a combined generating vector formed from two 
generating vectors a distance I apart; and 

b. datapoint loop means operating on each of the datapoints of 
said input signal for combining one partial possible message s E 
with a datapoint of said input signal and with an external 
portion of said combined generating vector thereby to update 
an element of said u vector defined by an internal portion of 
said combined generating vector. 

22. A method for synchronizing to pilot signals of a code division multiple 
access (CDMA) communication systems, the method comprising the 
steps of: 

a. per a set of possible PN loadings: 

i) loading a local PN generator with an initial PN loading; 
5 ii) looping on each of the datapoints of said received pilot 

signal, said step of looping including the steps of. 

a) combining one partial possible loading s £ with a 
datapoint of said received pilot signal and with a 
PN loading produced by said local PN generator 
2Q thereby to update a u vector and 

b) stepping said local PN generator to produce 
another PN loading; 
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Hi) performing a fast Hadamard transform on said u vector 
to determine the quality, in accordance with a metric, of 
each of said set of possible PN loadings; 
iv) selecting the PN loading having the best value for its 
metric; 

b. repeating steps i - iv for all sets of PN loadings; and 

c. selecting the PN loading with the best metric from among the 
PN loadings selected in step iv. 

23. A method for synchronizing to pilot signals of a code division multiple 
access (CDMA) communication systems which have frequency drift 
therein, the method comprising the steps of: 
a. per a set of possible PN loadings: 

i) looping over plurality of drift loop values, the step of 
looping including the steps of: 
g a) loading a local PN generator with a different initial 

PN loading per loop value; and 
b) generating an input signal from a received pilot 
signal which is insensitive to drift; 
ii) looping on each of the datapoints of said input signal, 
20 said second step of looping including the steps of: 

a) combining one partial possible loading s E with a 
datapoint of said input signal and with a 'PN 
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loading produced by said local PN generator 
thereby to update a u vector; and 
b) stepping said local PN generator to produce 
another PN loading; 

s Hi) performing a fast Hadamard transform on said u vector 

to determine the quality, in accordance with a metric, of 
each of said set of possible PN loadings; 
iv) selecting the PN loading having the best value for its 
metric; 

10 b. repeating steps i - iv for all sets of PN loadings; and 

c. selecting the PN loading with the best metric from among the 
PN loadings selected in step iv. 
24. A method for synchronizing to received encoded signals having 
frequency drift therein, the method comprising the steps of: 
15 a. per a set of possible PN loadings: 

i) looping over plurality of drift loop values, the step of 
looping including, the steps of: 

a) loading a local PN generator with a different initial 
PN loading per loop value; and 

b) generating an input signal from a received pilot 
signal which is insensitive to drift; 

ii) looping on each of the datapoints of said input signal; 
said second step of looping including the steps of: 
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a) combining one partial possible loading s £ with a 
datapoint of said input signal and with a PN 
loading produced by said local PN generator 
thereby to update a u vector, and 

b) stepping said local PN generator to produce 
another PN loading; 

Hi) performing a fast Hadamard transform on said u vector 
to determine the quality, in accordance with a metric, of 
each of said set of possible PN loadings; 

iv) selecting the PN loading having the best value for its 
metric; 

b. repeating steps i - iv far all sets of PN loadings; and 

c. selecting the PN loading with the best metric from among the 
PN loadings selected in step iv. 

25. A method for decoding signals having messages therein and being 
encoded with error correcting codes, the signals having frequency drift 
therein, the method comprising the steps of: 
a. per a set of possible messages: 

i) looping over plurality of drift loop values I, the step of 
looping including the steps of. 

a) creating a combined generating vector formed 
from two generating vectors a distance I apart; 
and 
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b) generating an input signal from a received pilot 
signal which is insensitive to drift; 
ii) looping on each of the datapoints of said input signal, 
said second step of looping including the step of. 

a) combining one partial possible message s E with a 
datapoint of said input signal and with an external 
portion of said combined generating vector 
thereby to update an element of said u vector 
defined by an internal portion of said combined 
generating vector 

iii) performing a fast Hadamard transform on said u vector 
to determine the quality, in accordance with a metric, of 
each of said set of possible messages; 

iv) selecting the message having the best value for its 
5 metric; 

b. repeating steps i - iv for all sets of messages; and 

c. selecting the message with the best metric from among the 
messages selected in step iv. 

26. A unit according to any of claims 1 - 21 substantially as described 
20 hereinabove. 

27. A unit according to any of claims 1 - 21 substantially as illustrated in 
any of the drawings. 
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28. A method according to any of claims 22 -25 substantially as 
described hereinabove 

29. A method according to any of claims 22 - 25 substantially as 
illustrated in any of the drawings. 
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